Compression modulus of macroscopic fiber bundles. 
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^ We study dense, disordered stacks of elastic macroscopic fibers. Tliese stacks often exhibit non- 

, linear elasticity, due to the coupling between the applied stress and the internal distribution of 

fiber contacts. We propose a theoretical model for the compression modulus of such systems, and 
illustrate our method by studying the conical shapes frequently observed at the extremities of ropes 
and other fiber structures. 
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Threadlike objects are the constituents of many natural and synthetic materials At the microscopic level, 

filaments such as actin, microtubules and other semi-flexible polymers control not only the elastic and viscous prop- 
erties of different biological architectures but also of polymer melts, surfactant solutions and gels Q. These are 
collective properties that can be understood from the individual characteristics of the filaments like the bending and 
: stretching elasticity or the friction coefficient j5|]. At such small scales, thermal disorder plays a key role that has been 
successfully described within the framework of statistical mechanics. For instance, the reptation theory for polymer 
melts and solutions , ^ or the description of rigid gel elasticity |^ successfully account for the collective macroscopic 
behavior, based on a few, microscopic parameters of the single polymers. Surprisingly, macroscopic threadlike systems 
have been much less studied, in spite of their abundance in the natural and synthetic realms: hair, wool, cotton and 
^ ' other natural fibers or Nylon strings, ropes and textiles or glass wool belong to this class of systems 2]. Here also 
' the elastic or friction constitutive parameters of the individual objects are well known. For instance, bending and 
stretching moduli! of human hair, wool. Nylon and steel fibers can easily be found in the literature. But theories that 
predict collective fiber properties such as the compressibility of a hair tress, or studies that explain what forces play 
, a role in a disentanglement action like combing or carding, are scarse [MIS' In this letter we argue that methods 
^ ' developed for fiber microscopic systems |0 are also useful at the macroscopic level, and draw an analogy between the 
O [ effect of thermal fluctuations and the effect of the intrinsic disorder carried by the spontaneous shapes of the fibers. 

We consider a two dimensional fiber bundle as illustrated in figure ^ Discussion in two dimensions allows for a 
simpler presentation of the methods, while retaining the main ingredients of the problem. Extension to the three 
dimensional case will be discussed later. The bundle is composed of a large number Af of elastic strands of length L, 
in a space of height d x J\f, where d is the average distance between two fibers. The shape of the fiber is described by 
the function Cni^) that measures heigth deviations from the line of average position at height d x n. This description 
in the Monge paramctrization is well adapted to the shapes of interest, with a weak gradient dQn{x)/dx ^ 1. In the 
, absence of interactions between fibers, each filament n has a spontaneous shape, described by the function (^o,n{a^)- It 
' is the disorder associated with the spontaneous shapes of the different fibers that leads to a non-trivial compression 
' behavior of the whole system. Under an external force, the spontaneous shape of the fiber Co,n(a;) will be transformed 
' into the actual shape C,n{x). Within the elastic limit, the deformation of a segment of length dx around the natural 
-j-J ' shape will cost an elastic energy d[/b, proportional to the bending modulus k and to the curvature differences, 
d[/b-^/2(C(:r)-C^',„(x))2dx. 

The bending modulus k is a constitutive property of the fiber, that depends on the shape of its section and on the 
Young modulus of the material 5j. For instance, for a cylindral filament of radius a and Young modulus E one has 
^ ■ K — Ea*/A. The total bending energy of the JV filaments can then be written as : 
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If one assumes that different fibers interact through a two body potential then the total interaction energy Ui 
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can be written as 
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Uint = V / dx dx' V [Crnix) + md - {Cn{x') + nd)] . 
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In order to bring the structure of the interaction energy to a tractable level, we follow a simple mean-field approxi- 
mation, similar to the Helfrich treatment of thermally activated fluid membranes and fibers [ij, jlj • We assume that 
forces between first-neighbors dominate the interaction energy, an exact assumption for excluded volume potentials 
and a good approximation for other short range forces such as those given by screened electrostatic potentials. Then 
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FIG. 1: Dense fiber bundle. The function (n{x) describes tlie sliape of tfie ntii fiber witli respect to a reference average level. 



the sum of the bending and interaction |(5J terms can be written as the effective energy 

Because of the weak gradients involved, inter-fiber interactions take place only at the same x coordinate. B is the 
strength of the mean- field quadratic potential, to be determined by the self-consistent prescription B = dxd'^ (e) /d(P 
where (e) = (i?eff) /{LAfd), is the average energy density. The equilibrium state of the fiber stack correponds to the 
minimum of the effective energy functional (|3J) with respect to the fiber configurations {Cn(2;)}, given the spontaneous 
shapes {Co,n{x)}, or more precisely, given the spontaneous curvatures {Co'n(2^)}- The spontaneous curvatures are to 
be treated as the random functions responsible for disorder in the system. Averages, denoted by (), will then be 
performed with the probability associated with the spontaneous curvature distribution. Before carrying forward the 
outlined procedure, it is of interest to review results for the thermal case, as this will help to build intuition for the 
key quantities controlling the stack structure. 

In thermal stacks or bundles of fibers, disorder is introduced by the random Brownian motion of the fibers, which 
would otherwise be flat, i.e. {Co,rt(a;) = 0}. The shape of a single fiber in solution is described by a persistent walk, a 
class of statistical walks where tangent correlations vanish exponentially over a persistence length £p, related to the 
fiber rigidity through ip — k/ {ksT). A persistent walk primarily moves along the forward direction for a short enough 
distance ^, but it performs also side excursions of size D. From the equipartition theorem one has IpD^ jl^ ~ 1. In 
a stack, lateral excursions are limited by near-neighbors setting D = d. Typically, there is one interaction between 
fibers over a collision distance £c — ip^^d'^^'^. The energy density of the thermal stack (bt), being of the order of fc^r 
per collision area, (ey) — ksT/ldic), one gets a compression modulus: 

1 /3 

^ m: ^ iknTr/\-y^d-'/\ (4) 

which holds for large enough rigidities, £p ^ d. 

In macroscopic fiber stacks, the disorder is quenched and introduced by the spontaneous disordered shapes (^o^nix). 
Lets consider a typical shape of amplitude Co a-nd wavelength Qq^ ■ Clearly, there are no interactions for distances 
larger than the shape amplitude, d ^ Co- For smaller distances the local gradients are of order of Co^O: leading to a 
collision length £c — d/(Co9o)- The energy density is thus given by (e) ~ nQq^id {die) — K.Qlo/d, independent of Ic- 
It follows that the compression modulus reads: 



(5) 
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It is worth at this stage to stress several striking differences and similarities between the thermal and the macroscopic 
cases. First, contrary to thermal stacks, where the compressibility is only governed by the persistence length £p and by 
the distance between fibers d, in macrosocopic fiber stacks the spacial scale associated with disorder inhomogeneities 
plays also an important role. As a matter of fact, the scale for compressibility is set by the combination kqq: stiffer 
fibers lead to harder macroscopic materials, quite at the opposite of thermal stacks where increasing fiber stiffness 
reduces fluctuations and softens the stack. Also, while the cost of conflncmcnt for thermally controlled fluctuations 
diverges as the conflnement distance vanishes, only a finite amount of energy is required to completely flatten a 
macroscopic fiber. But in both cases, the compression niodulii increase with the degree of disorder, measured by the 
temperature in the thermal system, and by the square amplitude of spontaneous shapes, Q for the quenched bundle. 
And finally, the power law dependence of the thermal modulus is larger than the macroscopic one: quenched disorder 
appears to lead to more "robust" materials as these are less susceptible to variations in stack density. We now return 
to a more detailed calculation of the compression modulus. 

It is convenient to perform the functional minimization of the effective energy Q in the space of the eigenfunctions 
{$g} for the biharmonic operator, {d^/dx'^ — q'^)^q{x) = 0. We choose also boundary conditions ensuring that 

there are no forces nor torques on the strands extremities: <I>^ (0) 
orthonormal set of eigenfunctions can be written as 



= and $^^^(0) = <^t'{L) = 0. The 
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4>q{x) = 



cosh(ap) — cos(ap) 
sinh(Q;p) — sin(Q!j,) 
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where the numerical coefficients ap are determined from the relation cos (ctp) cosh (ap) — 1. The solutions obey 



approximately ap 
eigenfunctions ©, 
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TT with p G N. We first develop the shapes (x) and Co,n {x) on the basis of the 
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where Q is a wavevector along the stack direction, then minimise with respect to the coefficients C,qQ and finally 
compute the average energy density: 
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where B{Q) — 2B{1 — cosQd). Equation ^ shows how the spatial inhomogeneities of spontaneous shapes determine 
the energy density of the stack: the harmonic, elastic nature of the energy penalty for shape distortion and interfiber 
interactions translates into the sole dependence of the energy density on the values of the second moment of the shape 
disorder distribution, weighted by the classical correlation kernel for bending elasticity. A practical consequence of 
equation ^ is that it allows to connect fiber geometric information, available from a simple statistical analysis on 
fiber shapes, with collective properties such as the compression modulus B. As an example, we now calculate the 
compression modulus for monomodal structures of uncorrelated fibers, for which the standard square deviation is 
written as (Co.ijQ^) = ^LQ6qq„. In the hmit where A/" » 1, the dimensionless compression modulus B — B/^hqq) 
obeys the self-consistent differential equation: 
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where Co now sets the natural distance units, d = d/Co- At short distances, the compression modulus is high, and 
equation IjlUI) leads to the asymptotic form discussed in 0, _B = nq'^Q/d'^. For distances larger than the natural 
amplitude for spontaneous shape disorder, Coi the fibers don't interact significantly, and the compression modulus 
vanishes exponentially B = KqQ{d/<^o)exp[—d/(^o]. A numerical solution of equation (|10|l interpolating between these 
two regimes is presented in figure 

We now illustrate our results by studying the shape of the cones that appear at the end of ropes and other fiber 
bundles. The cones are generated by the expansion of the fibers, under the repulsive field that we discussed previously. 
The repulsive field experienced by one fiber, U{d) = B{d)d^ /2, can be roughly approximated by a step function of 
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FIG. 2: Full line is the compression modulus i? of a fiber stack vs interfiber distance d. Dashed lines represent the asymptotic 
behaviours at short distances, Bd? ~ 1, for d <^ 1 and at large distances B ~ dexp[— d ], for d S> 1. Plots are for dimensionless 
quantities B = -B/(«go) and d — d/^o, with k the fiber rigidity, qo and the typical wavevector and amplitude of spontaneous 
shape disorder. 

amplitude and range Co- The average path of one fiber is a straight, horizontal line, as long as it remains 

inside the constrained zone, as depicted in figure When it escapes confinement, the fiber crosses an expanding 
zone where it still interacts with its neighbors. Under the repulsive forces, the fiber bends until it escapes the repulsive 
potential, for interfiber distances larger than the range of the potential Co- Bending of the fiber is caused by two 
localized forces of amplitude / = k^oCo/^ and opposite signs, applied by the two closest neighbors. After escaping the 
expanding interaction zone, the fiber follows a straight line, at an angle that depends on its distance from the center 
of the bundle. If one enumerates the fibers in the stack from n = 1 in the center to n = JV/2 for the external fiber, 
the following set of equations recursively determine the position [imUn] where a fiber escapes its neighbors potential 

{In-l - £n){3il_, - {(n-l - inf) = j (y„_l - (n - (11) 
i£n-l - tn)tl - j(yn ~ nd - {(„ - d)) , (12) 

where is measured from the border of the confinement zone and y„ from the center of the stack. The distance 
d is the interfiber separation in the confinement zone. Equations (|12|l and the boundary conditions yi = Co and 
^%'/2 ~ iy^f/2 ~ dJ\f/2)3K/ f determine the shape to the boundary yn{^n), shown also in figure (|2Jl. The extension of 
the expansion zone, ii, depends on the number of fibers and on the compression rate a — (Co ~ '^)/Coj but not on the 
fiber rigidity: £i ~ CoA/'^/'^a^/'^(qoCo)~''^'^- Similarly, the geometric shape of the cone, given by the slope of the last 
fiber Oj\f/2: does not depend on the fiber rigidity, it is only fixed by the geometric quantites associated with the fiber 
disorder, 9^f/2 ~ a^/^{qoCay^^^ ■ 




FIG. 3: Cones are formed at the unconstrained ends of ropes and other bundles of fibers, as shown on the left for a Nylon 
rope. On the right, the calculated shape in two dimensions. Here, only the average fiber path is represented for simplicity. The 
fibers only interact mutually in the expansion zone, to the left of the curved line, and follow a straight path otherwise. The size 
of the expansion zone and the angle of the cone are independent of fiber rigidity, they are fully determined by the geometric 
properties of fiber disorder. 

For high enough fiber density, the compression modulus of a three dimensional bundle can be computed in manner 
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similar to the two dimensional case. Indeed, if the density is so high as to provide a cage environment for each of 
the fibers, as depicted in figure (0)), the forces between fibers are controlled by excluded volume interactions between 
near-neighbors. In this case, we define a two dimensional fiber density in the bundle a = Mi:a? /A where A is the 
area of the bundle section. For a perfectly ordered hexagonal array one would have a = 27r\/3/3 a?r~^ , where r is the 
distance between fiber centers. The maximum close packing density is (Tmax — 0.91. The relevant fiuctuating distance 
rf = r — 2a is now a function of the fiber density a. In the hexagonal case, the cage would hold up to to distances 
d = 2a or densities as small as CTc — 7r-\/3/24 ~ 0.23. In this regime, the compression modulus is expected to vary as 
S3D ~ '''Q'oCo''"^- explained above, for the non-thermal case, only a finite energy is required to completely flatten 
a fiber. Thus, the confinement energy density is here bound to a maximum value < esD >'~-" '^ZoCo'^"^; contrary to 
the three dimensional thermal situation jlO|, where the energy density diverges as '-^ cP^^. For densities smaller than 
(Tc, the situation is unclear, even for the thermal case. Large lateral excursions of the fiber are now possible 13], and 
there is not yet, to our knowledge, a clear criterium to determine the maximum allowed size of the excursions. We 
hope to address this question in future work. 



FIG. 4: Cross section of a three dimensional bundle in the high density limit where each fiber is confined in a cage defined by 
its first nearest neighbors. 

As a summary, we have shown that stacks of macroscopic disordered fibers behave as a compressible material, 
the compression modulus being proportional to fiber rigidity, at the opposite of microscopic fiber stacks where it 
decreases with fiber rigidity. We have also shown how the compression modulus depends on stack density and on 
the geometric quantities characterizing spontaneous shape disorder. Our results allow to predict collective properties 
from the individual fiber characteristics and provide a new tool to interpret experiments in fiber systems. As an 
example of this we studied terminal bundle cones, showing that this class of shapes follows an universal behavior that 
depends only on fibers geometric features. Extension of our results to three dimensions is straightforward in the limit 
of large fiber density. The weak density limit poses an interesting challenge for fiber and polymer systems, both in 
the thermal and quenched cases. 
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